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We examined the rare decays B -> X s £+r and B s -> £ + 1 '7 in the minimal supersymmetric 
model with large tan f3. Taking into account the gluino-loop and neutralino-loop effects, we found 
that for a large tan/3 the neutral Higgs exchanging diagrams could enhance Br(B — > X s t + t~) by a 
factor of 5 and Br(B s — > t + t~7) by a couple of orders in some part of supersymmetric parameter 
space allowed by current experiments such as b — » sy, B — » K^£ + £~ and B s — > £ + £~. The 
forward-backward asymmetry and the distributions of differential branching ratios are also found to 
differ significantly from the standard model results. Such enhanced branching ratios reach the level 
of 10 -5 and thus might be observable in the new generation of B experiments. 
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I. INTRODUCTION 



Flavor-changing neutral-currents (FCNC) induced B-meson rare decays provide an ideal opportunity for extracting 
information about the fundamental parameters of the standard model (SM), testing the SM predictions at loop level 
and probing possible new physics. After the observation of the penguin-induced decay B — > X s "f and the corresponding 
exclusive channels such as B — > K*j M, rare B-decays have begun to play an important role in the phenomenology 
of particle physics. The latest measured decay ratio for B — > X s -f by CLEO and BELLE Q is in good agreement 
with the SM prediction, putting strong constraints on its various extensions and therefore stimulating the study of 
radiative rare B-meson decays with a new momentum. 

Among rare B-meson decays, B s — > £ + £~~{ (£ = e, /i, t) are of special interest due to their relative cleanliness and 
sensitivity to models beyond the SM Since in these processes a photon is emitted in additional to the lepton 

pair, no helicity suppression exists and "large" branching ratio is expected. Other interesting decay modes in this 
context are the inclusive transitions B — > X s £ + £~. Although these rare decays have not been observed, their detection 
is expected at the B-factories which are currently running. 

These decays have been studied in the SM |jj and recently, to reduce the theoretical uncertainties, the next-to-next- 
to leading order(NNLO) corrections were completed @. New physics effects in these decays have also been studied 
in some models, such as the minimal supersymmetric model (MSSM) |p[-pT[, the two Higgs doublet model(2HDM) 
^2|-|l5| and the technicolor model ||. 

It is noticeable that in the SM the matrix elements of B — > X s £ + £~ and B s — > £ + £~j are strongly suppressed by 
a factor mt/myy and the contributions from exchanging neutral Higgs boson can be safely neglected. In the MSSM 
fbsf , the situation is different, specially in the case of £ = r with large tan/3. In this model, as studied in p| , the 
contributions from exchanging neutral Higgs bosons are enhanced roughly by a factor tan 3 j3 and may no longer be 
negligible since currently a large tan /3 is favored both by LEP experiments |17| and by the supersymmetry (SUSY) 
explanation [ p"8| of the measurement for the muon anomalous magnetic moment p"9[ . 

We note that the most previous studies on the contributions from exchanging neutral Higgs bosons ]9|-|ll|| mainly 
focused on the charged-current loop effects; the effects of neutral-current loops (NCL), such as gluino-loop and 
neutralino-loop, have been considered only for B® s — > using mass insert approximation method in B. A 

detail general calculation for such NCL effects is necessary. The NCL can be induced via the flavor mixing of down- 
type squarks and might be important for the following reasons. Firstly, such flavor mixings of sfermions are almost 
unavoidable in supersymmetric models. In fact, in the framework of MSSM sfermions may have arbitrary flavor 
mixings in the soft breaking terms; while in some constrained MSSM, such as low-energy supergravity models, the 
flavor mixings at weak scale could be naturally generated through renormalization equation even the flavor diagonality 
is assumed at the Planck scale |20) . Secondly, the flavor mixings between the third and the second generation squarks 
are subject to no strong low-energy constraints like K° — K° mixing. Thirdly, the large tan/3 will give rise to large 
mass splitting between two mass eigenstates of sbottoms, making the lighter sbottom (b\) even lighter. 

In this article, we will present a complete calculation of MSSM effects in the decays B s — > £ + £~j and B — > X s £ + £~ 
(£ = e, fi, t), taking into account the contributions from the neutral Higgs exchange with NCL. We will evaluate 
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the effects on branching ratios, the forward-backward asymmetry as well as the distributions of differential branching 
ratios. In Section 0, we will give a brief description of the squark mixing in the MSSM. The detailed calculations 
for Wilson coefficients of scalar and pseudo-scalar operators are included in Section [II. The decays B — > X s £ + l~ 
and B s — > £ + £~j are calculated in Section |y|. Experimental constraints on the parameter space of th e MSSM are 
discussed in Section ^. Some numerical results are presented in Section [v|. Finally, in Section VII , we give our 
conclusion. 



II. SQUARK MIXINGS IN MSSM 

The squark mass terms arise from the scalar potential which contains the supersymmetric F-term and D-term as 
well as the soft SUSY breaking terms. These soft breaking terms may have arbitrary flavor mixings. As a result, the 
squark mass terms in flavor basis, i.e., U = (ul, cl, th, ur, cr, tn) and D = (<1l, sl, bL, Ar, sr, 6r), take the forms 
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1/3, mjj( D -} is the 3x3 diagonal mass matrix for up(down)-type quarks. M^, M% and 

:es 

Vijv\ is ratio of 



where ejj = 2/3, eu ■ 

M 2 ^ are soft-breaking mass terms for left-handed squark doublet Q, right-handed up and down squarks, respectively. 
A u (Ad) is the coefficient of the trilinear term H2QU (HiQD) in soft-breaking terms and tan/3 
the vacuum expectation values of the two Higgs doublets. The hermitian matrices M? ^ can be diagonalized by 
the unitary relations, which transfer the interaction eigenstates into the physical mass eigenstates. So in the general 
MSSM, without knowing the mechanism of SUSY breaking, squarks could have arbitrary flavor mixings. 

However, the flavor mixings in the first two generations are subject to strong phenomenological constraints, such 
as K° — K° mixing. So we only consider the flavor mixings between the second and third generations, i.e., between b 



and s. Further, like the analysis in |22|, we suppose the tree level Lagrangian is flavor diagonal and the flavor mixing 
is induced via loops. The dominant effects are from the logarithmic divergences caused by soft breaking terms. Such 
divergences must be subtracted using a soft counter-term at the SUSY breaking scale, such as Planck scale M p . Thus 
a large logarithm factor \n(M 2 /m 2 v ) w 80 remains after renormalization. In the approximation of neglecting the 
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In 



Here 9b is the left-right sbottom mixing angle, Ka are the CKM-matrix elements, and c;, — % 

u G tt & ' b J ' u 4-7T 2mw Sill - ^ n 

Note the mixing matrix T u between the left-handed scharm and stops takes the similar form as T D , and under the 
assumption that the flavor mixing between the second and third generation squarks is at least one order lower than 
the third left- and right-hand squarks mixing, the rotation matrix T U ' D are approximately unitarity. 
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III. CALCULATION OF SCALAR AND PSEUDO-SCALAR WILSON COEFFICIENTS 



In the MSSM, the short distance contribution to b — > s£ + £~ decay can be computed in the framework of the QCD 
corrected effective weak Hamiltonian, obtained by integrating out heavy particles, i.e., top quark, W^, Z bosons in 
the SM and the sparticles, 
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(3.1) 



where Oi, Qi are operators given in Eqs. (7J) |12|| , ( |7.2| ) |13| and Ci, CQ i are Wilson coefficients renormalized at the 
scale fi r R. 
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FIG. 1. The Feynman diagrams involving charged Higgs which give dominate contribution to Cq 1 



In the MSSM the additional contributions to operators in Eq. ( |3.l| ) can be characterized by the values of the 
coefficients Ci and Cg i at the perturbative scale tovv- For the processes we will study, it is only relevant with the 
effective Wilson coefficients C7.940 which have been computed in Ref. and Cq 1 2 of additional scalar and pseudo- 
scalar operators. In this section we will focus our attention on the calculation of the Wilson coefficients Cq 1 2 (mw) 
with the assumption that except for the third generation squarks all sfermions are degenerate and have a mass of ~ 1 
TeV. Besides the contribution from box diagrams, the neutral Higgs-bosons exchange diagrams include totally five 
classes of loops: (1) W boson and up-type quarks, (2) charged Higgs and up-type quarks, (3) charginos and up-type 
squarks, (4) neutralinos and down-type squarks, and (5) gluinos and down-type squarks. As pointed out in Section |, 
the last two classes of loops have not been calculated in the literatures. Now we take all of them into account and 
use the Feynman rules presented in |l6| . Since we are only interested in large tan/3 case, for simplicity we ignore less 
important terms and keep only the leading part contributions given in the following: 



Charged Higgs (Fig. ||) 



C^(mw) = -Cg{m w ) = (x H±t ) (3.2) 



' w 



x The most general Hamiltonian in low-energy supersymmetry also contains the operators O^, Q[ which are flipped chirality 
partners of Oi, Qi. However, they give negligible contributions and thus not considered in the final discussion of physical 
quantities E3[. 
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where Xij = mf / m 2 and the one-loop integral functions Pi are given in the appendix. 
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FIG. 2. The Feynman diagrams involving chargino which give dominate contribution to Cq 1 . 
Chargino (Fig. ||) 
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(3.3) 
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where a is the mixing angle of neutral components of the two Higgs doublets, r + = cos 2 am H 2 + sin 2 am h 2 , ro = 
sin2a(m^ 2 — m^ 2 ) At = m^/y2mw sin/3 and A& = mj/v^mw cos/3 are the Yukawa couplings of top and bottom 
quarks, respectively. E/, and iV are the matrices which diagonalise the chargino and neutralino mass matrices. 
m — mj = mg is defined as the common mass of the first two generation squarks and all sleptons which are assumed 
to be degenerate, N'a — 3 ta,n0wNji — Nj2, and 
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FIG. 3. The Feynman diagrams involving neutralino which give dominate contribution to Cq 1 

• Neutralino (Fig. ||) 
The formula for the contribution from the neutralino are given by 
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Gluino (Fig. 
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It is noticeable that the supersymmetric contributions to Cq x 2 have an overall enhancement factor tan 2 /3. Moreover, 
the gluino loop contributions have an additional enhancement factor , i?/"^* — ■ So sizable contributions from 
neutral Higgs penguin diagrams are expected for a sufficiently large tan (3. 



IV. B RARE DILEPTONIC DECAYS IN THE MSSM 



A. Inclusive decay B — > X s £ + £ 



Neglecting the strange quark mass, and with p standing for the momentum transfer, the effective Hamiltonian (3.1) 
leads to the following matrix element for the inclusive b — > s£ + £~ decay, 
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The Wilson coefficients can be evaluated from mw down to the lower scale of about m b by using the renormalization 
group equation. When evolving down to b quark scale, the operators Oj (j = 1 — 6), Q3 can mix with 0^, (i = 7, 9); 
however, they can be included in an "effective" ©7,9 because of their same structures contributing to the b — > s£ + £~ 
matrix element. As for long-distance contribution from the intermediate J/ty family we follow Ref. |25| and include 
the effect in "effective" C| . Expanding d in powers of a s , i.e., C — C° + f^C 1 , one obtains the leading order 
effective Wilson coefficients |I1^-|T4| 
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Here s = p 2 /m 2 is the scaled dilepton invariant mass square, , r\ = ot s (mw) 1 0- s {mb) , vector hi, dj, and are given 
in |1 and 
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At NLO level, the Wilson coefficient C\ ,e ^ can be found in 
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Function g(m 2 /r7i 2 , s) in Eq. (4.3) arises from the one-loop matrix elements of the four-quark operators, and 
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where 



Z(x,y) = 



2 4z 



(4.11) 



To estimate the long-distance contribution in the second term in brace of Eq. ([0^), we take the phenomenological 
parameter k as 2.3 [|25| in numerical calculations. 

The formula of invariant dilepton mass distribution has been derived in |l3j ] , which is given by 
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where r = m 2 /m 2 . To get rid of large uncertainties due to m\ and CKM elements in Eq. (4.12), we normalized the 
decay rate to the semileptonic decay rate 

T(B -> X c £v) 
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2Ax 4 lnx is the phase-space factor and k(x) is a sizable next-to-leading QCD 



Here f(x) = 1 - 8x 2 + 8x 6 - x 
correction to the semileptonic decay 26 

The angular information and the forward-backward (FB) asymmetry are also sensitive to the details of the new 
physics. Defining the forward-backward asymmetry as 



Afb(s) 



J d cos 9(d 2 T/dsd cos 9) — d cos 9 (d 2 T/dsd cost 



(4.15) 



d cos 9 (d 2 T/dsd cos 9) + J\ d cos 9 (d 2 T/dsd cos 9) 
where 9 is the angle between the momentum of B-meson and £ + in the center of mass frame of the dilepton, we obtain 
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B. Exclusive decay B s 
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Now let us turn to the rare radiative decay B s — » £ + £~-f. The exclusive decay can be obtained from the inclusive 
decay b — > si? + £~ 7 , and further, from 6 — >■ s£ + £~ . To achieve this, it is necessary to attach photon to any charged 
internal and external lines in the Feynman diagrams of b — > s£ + £~ . As pointed out in Ref . |!| , contributions coming 
from the attachment of photon to any charged internal line are strongly suppressed and we can neglect them safely. 
However, since the mass of ^-lepton is not much smaller than that of i? s -meson, in B s — > £ + £~j decay, the contributions 
of the diagrams with photon radiating from final leptons are comparable with those from initial quarks. When a photon 
is attached to the initial quark lines, the corresponding matrix element for the B — * £ + £~j decay can be written as 
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Here e u a nd g M are the four vector polarization and momentum of photo n, respectively; Gi, F{ are form factors p7j , f28[ . 
Eg. (4. 2l] ) can be obtained by multiplying p^ in both sides of Eq. ( 4.2C ) and using the equations of motion. From Eq. 
(4.21) one can see that the neutral scalars do not contribute to the matrix element Aii. 
When a photon is radiated from the final ^-leptons, the situation is different. Using the expressions 
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and the conservation of the vector current, one finds that only the operators Qi l2 and Og give contribution to this 
Bremsstrahlung part. The corresponding matrix is given by B 
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Here P Bs , f Ba are the momentum and the decay constant of the B s meson, p\, p 2 are momenta of the final £-leptons. 

Finally, the total matrix element for the B s — > £ + £~-f decay is obtained as a sum of the A4i and M. 2 . After 
summing over the spins of the £-leptons and polarization of the photon, we get the square of the matrix element as 
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It is obvious that the quantity \Ai\ 2 depends only on the scalar products of the momenta of the external particles. 
In this paper, we follow Ref. |3| and consider the photon in B s — > £ + £~"f as a hard photon and impose a cut on the 
photon energy E 1 0, which correspond to the radiated photon can be detected in the experiments. This cut requires 
E 1 > 5 m Bs /2 with S = 0.02. 

After integrating over the phase space and the lepton energy E\, we express the decay rate as 



2 When photon is soft, both processes of B s —> £ + £ 7 and B s — * £^ 
infrared singular terms in |7V/( 2 1 2 can be canceled exactly by the 0(a e 



should be considered together, and in this case, the 
virtual correction in _B, — > £ + £~ [HI. 
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where s = p 2 /m^ , f = m 2 /m 2 B , z = Z(r, s) takes the form given in Eq. (4.11) 



(4.28) 



V. EXPERIMENTAL CONSTRAINTS ON THE RELEVANT SUSY PARAMETERS 

Before scanning the relevant parameter space of MSSM, we assume (i) the masses of the particles are restricted to 
the sub-TeV regime and larger than the lower experimental bounds J2lj ; (ii) A-parameters are smaller than 3Mq [^9| ; 
(iii) Except for the third generation squarks, all sfermions are degenerate and have masses of ~ ITeV; (iv) The GUT 
mass relation Mi w Mij2 for gauginoes is used; (v) Only small flavor violation e,; < 0.1 is allowed. We also take into 
account the well-known large radiative corrections to neutral Higgs masses Q . With above assumptions and gluino 
mass fixed as its lower experimental bound 190 GeV [2lJ, the relevant parameter space is determined by nine input 
parameters: Mq, M^ r , M^^, fj,, tan/3, A t , A b , Mi and itia- In addition, we consider the following the experimental 
constraints in our scan: 

(1) The recently reported value of muon g — 2 [|l9| shows a 2.6 standard deviation from its SM prediction. The SUSY 
explanation of this deviation requires (i) fi > and (ii) large tan/3. In our calculation we assume 20 < tan/3 < 50 
and at least one of the charginos or neutralinos must be lighter than 500 GeV |Ts| ]. 

(2) Non-observation of any supersymmetric signals at CERN e + e~ collider LEP-II and the Fermilab Tevatron 
imposes lower bounds as 



m H ± > 78.6 GeV m h o > 88.3 GeV, m~± > 67.7 GeV, 
m^o > 42.0 GeV m~ ti > 86.4 GeV, m b [ > 75.0 GeV. 



(5.1) 



(3) The latest measurement of the inclusive branching ratio from CLEO and BELLE Q gives world average value 

2.44 x 10~ 4 < Br(B -> X sl ) < 4.02 x 10~ 4 (95% C. L.), (5.2) 

which is specially useful to constrain extensions of the SM. Previous studies used leading-order (LO) SM result to limit 
the MSSM parameter space. Instead of LO calculation p^| , pl| , the branching ratio of B(B — > X s j) has been estimated 
at next-leading-order (NLO) level in the SM with about 22.5% increase of its central value and theoretical error 
less than 10%. In this paper, we use the SM result computed at NLO level whereas the additional supersymmetric 
contributions at LO level. We do not use the available NLO matching conditions for the supersymmetric particles 
since they are computed under the specific assumptions about the sparticle spectrum, not necessarily satisfied in the 
criteria, and moreover, they are not valid for large values of tan/3 f|. Since the Wilson coefficient C% accounts for 
only 3% of the standard model b —> amplitude, it is therefore not expected to be significantly more important in 



the MSSM. In this case, Eq. ( |5.2| ) implies the constraints on the ratio R7 — G^'^jgg^j /C^' e ^ ' 



SM 

0.83 < R 7 < 1.13, or - 1.24 < R 7 < -0.94. (5.3) 

(4) The stringent bounds on the magnitude of the short distance coefficients come from the Collider Detector at 
Fermilab (CDF) |H| 

Br{B+ -> K+e + £-) < (60, 5.2) x 10~ 6 (90%C.L.), 
Br(B° -> K*°e+r) < (290,4.0) x 10~ 6 (90%C.L.), 

Br(B s — > £+£- ) < (5.4,2.0) x 10" 6 (95%C.L.) (5.4) 

for £ = e,fi. 
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VI. NUMERICAL RESULTS AND DISCUSSIONS 



We perform a complete scan over the nine-dimensional parameter space of the MSSM. For reference, we present 
our SM predictions 

Br(B -> X s £ + £~) = (11.8, 9.87, 2.94) x 1CT 6 , (6.1) 
Br(B s -> = (2.50, 2.62, 5.49) x 1CT 8 (6.2) 

for I = e,fi,T. These values are obtained by taking the QCD coupling constant a s (mb) — 0.218, the masses, decay 
widths and branching ratios of J/\& family in Ref. pi] , the normalized factor Br(B — > X c lv) = 10.2% and the 
form-factors [^8| 

2 1 GeV 2 3.74 GeF 

Gl(P } " (l-p 2 /5.6 2 ) 2 ' G2(P } " (l-p 2 /40.5) 2 ' 
„ , 2 , 0.8 GeF „ . 2N 0.68 GeV 

Fl(p)= (l-pV6-5 2 ) 2 ' F2b)= (l- P W ( } 

as well as the fixed input parameters pl[ listed in Table I in the numerical calculations. 



Table I. The value of the input parameters used in the numerical calculations (mass and decay constant in unit GeV). 



m t 


m c 


m b 


m T 






176 


1.4 


4.8 


1.78 


5.26 


80.45 


Ib s 


\K tb K* s \ 


\K tb K* ts /K cb \ 2 


a" 1 

u em 


r(5 s ) 


sin^ 0w 


0.21 


0.045 


0.95 


137 


1.64 x 10~ 12 s 


0.233 
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The results shown in Fig. ^ indicates that in some part of parameter space, the enhancement factors R of the branching 
ratios can be 5 for the inclusive decay B — > X s t + t~ and a couple of orders of magnitude over the standard model 
predictions for the exclusive decay B s — > t + t~7. This is quite different from the previous studies Jo]. As illustrative 
examples, the branching ratios for ditau final state as functions of the gluino mass are plotted in Fig. ||. The effects 
drop with the increase of gluino mass, showing the decoupling property of the MSSM. Table 2 shows the case in which 
both enhancement factors reach their maximum values and the corresponding SUSY parameters. A comparison of 
the enhancement factors with and without the NCL effects is also presented. 



Table II. The maximum enhancement factors and their corresponding SUSY parameters (mass in unit GeV). 



Rmax(B -> X s tt) 


Rmax{B s -> TTj) 


With NCL 
4.34 


without NCL 
0.99 


With NCL 
327.0 


without NCL 
0.75 


tan /3 = 40 m A = 453 = 143 


m Si = 92 



We have the following comments on the results: 

(1) The constraint on B s — > £ + £~ is of special useful to li mit th e effects of the scalars and pseudo-scalars. However, 
since no helicity suppression compared with Cio (see Eq. 4.23), the contribution of Cq 12 in B — > X s t + t~ and 
B s — > £ + £~j decay is still dominate for large tan/3. 

(2) As helicity suppression, the effects of Cq 1 2 in B — > X s t + t~ are very different from those with other dilcpton 
final states (see Eq. 4.13). The smaller mass of lepton in the final states, the less effect of gluino and the neutralinos. 
In the inclusive decays with lighter dilepton final states, the dominate contribution comes from the interaction term 
2Re(C? ff Cg ff ) with opposite sign to its SM value. For large tan/3 region, relative large m T make the neutral Higgs 
with NCL contribute sizably to the inclusive decay. 



'o 
o 

33 

2 

■c 
o 
c 

2 

CQ 



100 



400 



700 



1000 



Mass of gluino (GeV) 



FIG. 5. The branching ratios vs gluino mass with tan/3 = 40, tua = 271 GeV, mj 1 = 320 GeV, mg = 265 GeV. The solid 
and dashed lines correspond to the inclusive and exclusive decays with £ = r, respectively. 



Some numerical examples in the region of Rj > are presented in Figs. 0-0. Fig. |] and Fig. show the branching 
ratios of B — > X s t + t~ and B s — * t + t~j, respectively. For the specified parameter values in the figures, one sees 
that the contributions from gluino and neutralinos are dominant. When the mass of the CP-odd Higgs is less than 
700 GeV or tan/3 > 40, the branching ratio of B — > X s t + t~ can be enhanced by a factor 2, and for B s — > t+t - 7, by 
one order over the SM results. 
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^. 6 
5 



20 
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tan$ 



FIG. 6. Branching ratio of B 



X s t + t versus itia in (a) and versus tan/3 in (b). The dotted line stands 



for the SM prediction, whereas the solid (dashed) one is the MSSM prediction with (without) the contributions of 
the gluino and neutralinos. The parameter space in (a) is specified as m,- = (179,304) GeV, mi = (257,700) GeV, 
m x ± = (94,417) GeV,m^o = (81,109,200,417) GeV and tan = 20; In (b), m A = 899 GeV, the masses of the sparticles are 
dependent on the value of tan j3. 




100 400 700 1000 20 30 40 50 

m A (GeV) fanp 

FIG. 7. The same as Fig. |, but for B s -> r+T^-y. 

The differential branching ratios of B — > X s t + t~ and B s — > t + t~7 versus the scaled invariant dilepton mass 
squared are plotted in Fig. H, while the dependence of the FB asymmetry for B — > X s t + t~ on the scaled invariant 
dilepton mass squared is shown in Fig. || The figures show significant differences between the SM and the MSSM 
predictions, especially in large invariant dilepton mass region. 
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0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 

A 

s s 

FIG. 8. The same as Fig. ^a, but for the differential branching ratios of B — > X s t + t~ and B a —* t + t~j versus the scaled 
invariant dilepton mass squared with tua = 305 GeV. 

We should point out that since some common contributions appearing in both the numerator and the denominator 
cancel out to some extent, the FB asymmetry is a sensitive, relatively model-independent probe of these models. 
We stress that all these distributions would be useful for fitting the future experimental results in the framework of 
the MSSM, especially when some deviations from the SM predictions are discovered in future experiments. Different 
models, such as the MSSM and 2HDM's, may all predict some enhancements, but they may give different behaviors 
for some distributions. 




■0.5 - 



0.5 0.6 0.7 0.8 0.9 1 

s 

FIG. 9. The same as Fig. ^|a, but for the forward-backward asymmetry of B — > X s t + t~ versus the scaled invariant dilepton 
mass squared s with tua = 305 GeV. 

Although decay modes with the ditau final states are experimentally difficult compared to their di-muonic counter- 
parts, the considered decay modes are more sensitive to new physics. More theoretical studies, such as higher order 
effects which might also lead to important modifications as in the case of B — * X s ^ |52|, and experimental efforts for 
the decays considered are valuable. 
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VII. CONCLUSION 



In this work we performed a complete one-loop calculation of the inclusive decay B — > X s £ + £ and exclusive decay 
B s — > ^+^7 in the MSSM. Various experimental constraints on the relevant SUSY parameters, such as B — > A s 7, 
£? — > K^i + £~ , B s — ► and the latest g^ — 2 experimental measurement, were considered to constrain the 

parameter space of the MSSM. Our results showed that the contributions from the gluino and neutralino loops, which 
were neglected in previous studies, might be quite important or even dominant in some part of parameter space. 
These supersymmetric contributions could significantly enhance the branching ratios over the SM predictions. Also, 
with these contributions the distributions of the forward-backward asymmetry of B — > X s t + t~ and some other 
distributions could differ significantly from their SM predictions. Up to now the currently running B-factories such 
as BaBar at SLAC ||| and BELLE at KEKB @ have collected about 3.2 x 10 7 BB pair with a luminosity of 
(3 - 4) x 10 33 cm- 2 " [HQ. And the BaBar will take 10 8 BB pairs in three years while the BELLE 10 7 BB 
pairs each year. Since the branching ratios of the inclusive decay B — > X s t + t~ and exclusive decay B s — > t + t _ 7 can 
be enhanced by the SUSY effects with large tan/3 to reach the level of 1CU 5 , these ditau decays might be observable 
at the B factories. 
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APPENDIX 



In the Appendix we present the operator basis Oi, Qt and one-loop integral functions for b — > s£ + £ 
The operator basis 0{ is the same as the one used for the b — ► s£ + £~ in the SM, 

01 - {s a ^Lc ){c^Lb a ), 

2 = (a a 'fLc a )(cp'fLbf,), 



3 . 5 = (s a rLc a ) [ ^q P r{L,R)b p 



O4.6 = 



(s a ^Lc a ) I £\ 

V 9 

\ q 



0? = -r^s a a fJ,v (m b R + rrisL^aF^, 

l07T z 

°8 = -Shs a ^{m b R + m s L)t af} bpG%, 
° 9 = TT^(^Lb a )(£^£), 

1D7T 

do = T^(s a ^Lb a )(£^ l5 £), (7.1) 

where the chiral structure is specified by the projectors L, R = (1 T 75)/2, while a and [3 are color indices. and 
G" v denote the QED and QCD field strength tensors, respectively. t a p are the color triplet generators, g and g s stand 
for the electromagnetic and strong coupling constants. 

Operators Qi come from exchanging the neutral Higgs bosons in MSSM and are defined by 13 

2i = ^(s*Rb a )(££), 

lon z 

Q2 = -f^(s Q i?M(^), 

l07T Z 

„2 



15 



Q 5 ,6 
Q 7 ,8 



16tt 2 
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{s a Rbp) ^2sp(R 7 L)b a , 
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16 



16tt 

\ 9 / 

The one-loop integral functions which appear within the MSSM matching conditions are given by 

l 



(7.2) 
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(7.3) 



[1] R. Ammar, et al, CLEO Collaboration, Phys. Rev. Lett. 71 (1993) 674; M. S. Alam est al, CLEO Collaboration, 

Phys. Rev. Lett. 74(1995) 2885. 
[2] S. Glem, et al, CLEO Collaboration, |hep-ex/0104045| K. Abe, et al, BELLE Collaboration Phys. Lett. B511, 151, (2001). 
[3] T. Goto, Y. Okada, Y. Shinizu and M. Tanada, Phys. Rev. D55, 4273 (1997); T. M. Aliev, A. Ozpineci and M. Savci, 

Phys. Rev. D55, 7059 (1997); T. M. Aliev, N. K. Pak and M. Savci, Phys. Lett. B424, 175 (1998). 
[4] E. O. Iltan and G. Turan, Phys. Rev. D 61, 034010 (2000). 

[5] Z. Xiong and J. M. Yang, Nucl. Phys. B602 289 (2001); Y. Su, Phys. Rev. D 56, 335 (1997). 

[6] M. Misiak, Nucl. Phys. B393, 23 (1993); A. F. Falk, M. Luke and M. Savage, Phys. Rev. D49, 3367 (1994); 
Z. Ligeti and M. B. Wise, Phys. Rev. D53, 4937 (1996); F. Kriiger and L. M. Sehgal, Phys. Lett. B380, 199 (1996); 
A. Ali, G. Hiller, L. T. Handoso and T. Morozumi, Phys. Rev. D55, 4105 (1997); J-W. Chen, G. Rupak and J. Savage, 
Phys. Lett. B410, 285 (1997); G. Buchalla, G. Isidori and S. J. Rey, Nucl. Phys. B511, 594 (1998); G. Buchalla, G. Isidori, 
Nucl. Phys. B525, 333 (1998). 

[7] C. Bobeth, M. Misiak and J. Urban, Nucl. Phys. B574, 291 (2000); H. H. Asatryan, H. M. Asatrian, C. Greub and 

M. Walker, Phys. Lett. B507, 162 (2001). 
[8] P. H. Chankowski and L. Slawianowska, Phys. Rev. D 63, 054012 (2001); 

[9] P. Cho, M. Misiak and D. Wyler, Phys. Rev. D54, 3329 (1997); Y. Grossman, Z. Ligeti, and E. Nardi, Phys. Rev. D55, 
2768 (1997); J. L. Hewett and J. D. Wells, Phys. Rev. D55, 5549 (1997). 
[10] S. Bertolini, F. Borzumati, A. Masiero, and G. Ridolfi, Nucl. Phys. B353, 591 (1991); A. J. Buras and M. Miinz, 
Phys. Rev. D 52, 186 (1995); M. Ciuchini, G. Degrassi, P. Gambino, and G. F. Giudice, Nucl. Phys. B527, 21 (1998); 
B534,3 (1998). 

[11] C. S. Huang, W. Liao, Q. S. Yan and S.-H. Zhu, Phys. Rev. D59, 011701 (1999); C. S. Huang and and S.-H. Zhu, 
Phys. Rev. D61, 015011 (2000); C. S. Huang W. Liao, Q. S Yan and S.-H. Zhu, Phys. Rev. D63, 114021 (2001); 
C. S. Huang, W. Liao, Q. S. Yan and S.-H. Zhu, hep-ph/0110147t S. R. Choudury, N. Gaur, Phys. Lett. B451, 86 (1999). 

[12] B. Grinstein, R. Springer, and M. B. Wise, Phys. Lett. B202, 138 (1988); Nucl. Phys. B339, 269 (1990). 

[13] Y. B. Dai, C. S. Huang and H. W. Huang, Phys. Lett. B390, 257 (1997); 

[14] J. L. Hewett, Phys. Rev. D53, 4964 (1996). 



[15] H. E. Logan and U. Nier ste, Nucl. Phys. B 586, 39 (2000); C. Bobeth, T. Ew erth, F. Kriiger and J. Urban, |hep-ph/0104284| ; 

G. Erkol and G. Turan, ^ep-ph/011001~7| ; G. K. Yeghiyan, ^iep-ph/010815l| . 
[16] H. E. Haber and G. L. Kane, Phys. Reports. 117, 75 (1985). 
[17] See, e. g., the L3 collaboration, Phys. Lett. B503, 21 (2001) 

[18] L. Everett, C. L. Kane, S. Rigolim and L.-T. Wang, Phy. Rev. Lett. 86, 3484 (2001), J. L. Feng and K. T. Matchev, 
Phy. Rev. Lett.86, 3480 (2001). 



16 



[19] H. N Brown, et al., Mu g-2 Collaboration, Phys. Rev. Lett. 86, 2227 (2001). 
[20] See, e.g., M. J. Duncan, Nucl. Phys. B 221, 285 (1983). 
[21] Particle Physics Group. Eur. Phys. J. C, 15, 274 (2000). 

[22] K. I. Hikasa and M. Kobayski, Phys. Rev. D36, 724 (1987); J. M. Yang and C. S. Li, Phys. Rev. D49, 3412 (1994). 
[23] E. Lunghi, A. Masiero, I. Scimemi, L. Silvestrini, Nucl. Phys. B568, 120 (2000). 

[24] K. Chetykin, M. Misiak, and M. Miinz, Phys. Lett. B 400, 206 (1997); C. Greub and T. Hurth, Phys. Rev. D 56, 2934 

(1997). F. M. Borzumati and C. Greub, Phys. Rev. D58, 074004 (1998). 
[25] A. Ali and C. Greub, Z. Phys. C 49, 431 (1991); Phys. Lett. B 259, 182 (1991); 361, 146 (1995); A. Ali, B Ball, 

L. T. Handoko and G. Hiller, Phys. ReV. D61, 074024 (2000). 
[26] N. Cabibbo and L. Maiani, Phys. Lett. B79, 109 (1978); Y. Nir, Phys. Lett. B221, 184 (1989). 

[27] G. Buchalla and A. J. Buras, Nucl. Phys. B400, 225 (1993); D. Du, C. Liu and D. Zhang, Phys. Lett. B 317, 179 (1993). 

[28] G. Eilam, I. Halperin and R. R. Mendel, Phys. Lett. B361, 137 (1995). 

[29] H. E. Haber, Lectures Given at Theoretical Advanced Study Insstitute, 1992. 

[30] A. Djouadi, J. Kalinowski and M. Spira, Comput. Phys. Commun.108, 56 (1996). 

[31] M. Ciuchini, E. Franco, G. Martinelli, L. Reina, and L. Siliverstrini, Phys. Lett. B 316, 127 (1993); M. Ciuchini, E. Franco, 
G. Martinelli, L. Reina Phys. Lett. B 301, 263 (1993), Nucl. Phys. B 415, 403 (1994); M. Ciuchini, E. Franco, L. Reina, 
and L. Siliverstrini, Phys. Lett. B 421, 41 (1994); 

[32] G. Degrassi et al, JHEP0012, 009 (200 0); M. Carena et al., Phys. Lett. B99, 14 (2001). 

[33] BaBar Collaboration, B. Aubert et al, [hep-ex/0105044| ; Phys. Rev. Lett. 86, 2515 (2001). 

[34] Belle Collaboration, A. Abashian et al, Phys. Rev. Lett. 87, 091801 (2001). 



17 



